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ABSTRACT 

Recent X-ray observations reveal growing evidence for heating by active galactic nuclei (AGN) in 
clusters and groups of galaxies. AGN outflows play a crucial role in explaining the riddle of cooling flows 
and the entropy problem in clusters. Here we study the effect of AGN on the intra-cluster medium in a 
cosmological simulation using the adaptive mesh refinement FLASH code. We pay particular attention 
to the effects of conductivity and viscosity on the dissipation of weak shocks generated by the AGN 
activity in a realistic galaxy cluster. Our 3D simulations demonstrate that both viscous and conductive 
dissipation play an important role in distributing the mechanical energy injected by the AGN, offsetting 
radiative cooling and injecting entropy to the gas. These processes are important even when the transport 
coefficients are at a level of 10% of the Spitzer value. Provided that both conductivity and viscosity are 
suppressed by a comparable amount, conductive dissipation is likely to dominate over viscous dissipation. 
Nevertheless, viscous effects may still affect the dynamics of the gas and contribute a significant amount 
of dissipation compared to radiative cooling. We also present synthetic Chandra observations. We show 
that the simulated buoyant bubbles inflated by the AGN, and weak shocks associated with them, are 
detectable with the Chandra observatory. 

Subject headings: galaxies: active - galaxies: clusters: cooling flows - X-rays: galaxies 



1. INTRODUCTION 

Clusters of galaxies are excellent laboratories for study- 
ing the interaction between outflows from active galactic 
nuclei with diffuse baryons. Recent observational evidence 
demonstrates that the lives of AGN and galaxy clusters 
in which they reside are closely intertwined. In particu- 
lar, there is mounting observational evidence that AGN 
may provide a vital clue in explaining the cooling flow 
problem. In the absence of non-gravitational heating, the 
cluster cores should cool and accrete gas at rates of hun- 
dreds and more solar masses per year. However, this is in 
conflict with observational evidence which often indicates 
mass deposition rates of only ~ 10 solar masses per year. 
Moreover, the gas temperatures in cluster centers are typ- 
ically maintained above ~ 2 keV. As many cooling flow 
clusters are known to harbor active radio sources (Burns 
1990, Eilek 2004) and the enthalpy of cavities inflated 
by radio galaxies scales with the cooling flow luminosity 
(Birzan et al. 2004), AGN may serve as a viable heating 
source to prevent the gas from cooling and accreting at ex- 
cessive rates. Direct evidence for AGN heating has indeed 
been found in recent studies. Observations of the Perseus 
cluster (Fabian et al. 2003a, b) and the Virgo cluster (For- 
man et al. 2004) reveal sound waves and weak shocks 



in the intra-cluster medium (ICM). Even more recently, 
Nulsen et al. (2004) found evidence for shock heating 
in Hydra A, McNamara et al. (2004) in MS0735.6+7421 
and Sanderson, Finoguenov & Mohr (2004) in Abell 478. 
These results strongly suggest that AGN outflows can heat 
the ICM in a spatially distributed fashion, which may help 
to maintain ICM stability against radiative cooling. 

The evidence for non-gravitational heating in clusters 
has also been observed in a statistical sense in cluster scal- 
ing relations. These relations show departures from the 
self-similar scalings: In the absence of non-gravitational 
heating and radiative cooling, one would expect the en- 
tropy to scale proportionally with the mean cluster tem- 
perature. However, observations by Ponman, Sanderson & 
Finoguenov (2003), Pratt & Arnaud (2005) and Piffaretti 
et al. (2005) indicate a scaling of entropy roughly accord- 
ing to T^/^. Moreover, they reveal a systematic excess of 
entropy at large radii in low-mass clusters (e.g., Ponman, 
Sanderson & Finoguenov 2003). More observational evi- 
dence for cluster heating conies from a number of massive 
clusters that have been shown (McCarthy et al. 2004) to 
host very large cores which require an enormous amount 
of heating in order to explain their properties. In a recent 
study, Croston et al. (2004) separated a sample of groups 
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into radio quiet and radio loud objects. They demon- 
strated that radio loud groups deviate from self-similatity. 
This result demonstrates that AGN play a crucial role in 
heating the ICM and may offer the solution to the entropy 
excess problem. 

Based on observations of the Perseus cluster it first has 
been suggested by Fabian et al. (2003a, b) that viscosity 
may play an important role in dissipating energy injected 
by the central AGN. This idea has been tested in nu- 
merical simulations by Ruszkowski et al. (2004a,b) and 
Reynolds et al. (2005). Apart from viscosity, thermal 
conduction has been considered by a number of authors as 
means of transferring energy from hot outer layers of clus- 
ters towards the cool cluster cores (e.g., Fabian, Voigt, & 
Morris 2002, Voigt & Fabian 2004, Narayan & Medvedev 

2001, Briiggen 2003). 

The above observational results motivate our efforts to 
study the interaction between AGN and the ICM in de- 
tail. Numerical simulations of hot, underdense bubbles in 
clusters of galaxies have been performed by a number of 
authors (e.g. Churazov et al. 2001, Quilis et al. 2001, 
Briiggen et al. 2002a,b, Reynolds, Heinz & Begelman 

2002, Ruszkowski, Briiggen & Begelman 2004a,b, Dalla 
Vecchia et al. 2004). Common to these simulations is that 
they use a spherically symmetric, analytical profile for the 
ICM. Here, in order to get a more realistic representa- 
tion of the dissipative motions in the ICM, we perform 
hydrodynamical simulations of AGN heating in a cluster 
that has been extracted from a cosmological simulation. 
Our simulations include a dynamic dark matter that is 
represented by coUisionless particles and more sophisti- 
cated gas physics. The primary objectives of the current 
paper are to study the effect of viscosity and conduction 
on dissipation of AGN-induced gas motions and to assess 
the detectability of the signatures of AGN heating in syn- 
thetic Chandra observations. 

The paper is organized as follows. In Section 2 we dis- 
cuss the initial conditions and give the basic informations 
about the setup of the simulations. In Section 3, we discuss 
the implementation of heating by conduction and viscosity. 
Section 4 presents the method used to simulate Chandra 
observations. The results are discussed in Section 5 and 
summary in Section 6. 

2. INITIAL CONDITIONS AND THE SIMULATION SETUP 

The initial conitions for our simulations were computed 
with the SPH code GADGET in standard ACDM cosmol- 
ogy {flA — 0.7, Jim = 0.3, h = 0.7). They are based on 
a re-run of the S2 cluster in Springel et al. (2001) which 
we chose because of all simulated clusters available to us 
its properties were closest to those of the Perseus cluster 
(in terms of temperature profile, mass and baryon frac- 
tion). At redshift z = this cluster has a mass of 7 • 10^** 
M0 and a central temperature of ~ 6 keV. The initial 
density (right panel) and temperature (left panel) slices 
through the cluster center are shown in Figure 1. In Fig- 
ure 2 we show the initial entropy distribution. One can 
clearly see that this cluster is quite dynamic and shows a 
lot of substructure. It can also be seen in the temperature 



map that the central part of the cluster shows a cool core 
characteristic of cooling flow clusters. 

The SPH simulation of the cluster includes radiative 
cooling and star formation. The output of the SPH simu- 
lation serves as initial model for our adaptive mesh refine- 
ment (AMR) simulation. We use the FLASH code which is 
a modular block-structured AMR code, parallelised using 
the Message Passing Interface (MPI) library. It solves the 
Riemann problem on a Cartesian grid using the Piecewise- 
Parabolic Method (PPM) and, in addition, includes parti- 
cles that represent the coUisionless dark matter. The par- 
ticles are advanced using a cosmological variable-timestep 
leapfrog-method. Our simulations included 714346 col- 
lisionless particles that represent stars and dark matter. 
Gravity is computed by solving Poisson's equation with a 
multigrid method using isolated boundary conditions. For 
the relatively short physical time of the bubble simula- 
tion, radiative cooling and star formation are neglected, 
even though they were included in the cosmological SPH 
simulation with which the cluster was produced. The cen- 
tral cooling time in our cluster is about 700 Myrs which is 
much longer than the 140 Myrs of our simulation. So this 
is probably a fair assumption. The computational domain 
of our AMR simulation is a cubic box of side L = 2h~^ 
Mpc. For our grid, we chose a block size of 16^ zones and 
used outflow boundary conditions. The minimal level of 
refinement was set to 3 which means that the minimal grid 
contains 16 ■ 2*^"''^^' = 64'^ zones. The maximum level of re- 
finement was 7, which corresponds to an effective grid size 
of [16 • 2^^"^^]^ = [1024]^ zones or an effective resolution of 
1.96h~^ kpc. This resolution was chosen in order to cap- 
ture the relevant physical scales. Comparison with a run 
at a maximum refinement level of 8 (which we could afford 
only for a much shorter time) showed relatively little differ- 
ence in density and temperature and thus we think that we 
are close to convergence. Higher resolution tends to yield 
steeper waves and thus higher dissipation. Therefore, one 
may be confident that this simulation does not overesti- 
mate the amount of dissipation in the core. The code was 
run on 64 processors of the IBM p690 at the John-von- 
Neumann Institut for Computing in Jiilich, Germany and 
at the National Center for Supercomputing Applications 
at the Univeristy of Illinois on an identical machine. 

3. HEATING 

3.1. The energy and momentum equations 

The evolution of internal energy is followed by solving 
the energy equation 

de ^ . d f dT\ 

'di = -P^ + ''-'-+d^A"d^J' ^ ^ 

where k = 5.0 x 10'"^ {\n A/S7)-^T^/^ fc (Spitzer 1962) is 
the conductivity coefficient and fc is the conductivity sup- 
pression factor. The dissipation of mechanical energy due 
to viscosity, per unit mass of the fluid, is given by (Batch- 
elor 1967, Shu 1992, Landau & Lifshitz 1997) 

_ 2/z / 1 2\ 

^visc — — y^ij^ij ~ 3 ) ' ' •' 

where A = en and 





Fig. 2. — The initial entropy distribution in a slice through the center of the cluster. The contours are logarithmic and the colorbar indicates 
the logarithm in keV cm^. 
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dvi 
dxi 



(3) 



and where fi is the dynamical coefficient of viscos- 
ity. We use the standard Spitzer viscosity for an un- 

magnctizcd plasma (Braginskii 1958), for which /i = 

6.0 X 10-i^(lnA/37)-iT5/V« g cmT^ s^^, where is 
the viscosity suppression factor. We assume that both 
and /c are equal to 0.1. It is not certain that the sup- 
pression factors of both transport processes should be the 
same. Whether they arc the same may depend on the scale 
magnetic fluctuations extend to. This scale may be much 
larger than the gyroradii of electrons and ions (in which 
case suppression factors could be comparable) or it could 
be comparable to the ion gyroradius (E. Zweibcl, private 
communication). Nevertheless, we chose the same sup- 
pression factors for the sake of simplicity and the lack of 
a better choice. The magnitude of the suppression factor 
is motivated by various theoretical arguments (e.g., given 
in Narayan & Medvedev (2001) or by MHD simulations 
by Maron et al. 2004). However, we note that the precise 
value of the suppression factor is highly uncertain and, de- 
pending on the nature of magnetic turbulence, may even 
exceed the Spitzer value (Cho et al. 2003) or be supressed 
well below it. The supression factors may also vary as a 
function of distance from the cluster center. An observa- 
tional constraint on conductivity (/c < 0.15) comes from 
the fact that clusters would have cooled substantially over 
a Hubble time if conduction had not beed supressed (Loeb 
2002). 

As conditions inside the buoyantly rising bubbles are 
very uncertain and because we want to focus on energy 
dissipation in the ambient ICM, we assume that dissipa- 
tion occurs only in regions outside the bubbles. To this 
end we impose a condition that switches on viscous and 
conductive effects provided that the fraction of the in- 
jected gas in a given cell is smaller than 10^^. 

Velocity diffusion was simulated by solving the momen- 
tum equation 



d{pvi) d ^ dP , d'^ik 

o, + -^{pVkVi) + ^ =P9i+ ^ : 

at oxk oxi oxk 



where 



d 



2m ( eife - -A6ik 



(4) 



(5) 



and where all other symbols have their usual meaning. 

3.2. Energy injection by AGN 

The AGN is assumed to sit in the centre of the clus- 
ter. The intermittency period of the AGN, i.e. the period 
between two subsequent bursts, was set to 3 • 10^ years, 
within which the AGN was active only for a period of 5-10^ 
years. This choice was made in order to produce a number 
of strong waves in quick succession in order to gauge the 
dissipated energy within a feasible simulation time. The 
time-averaged luminosity of the AGN was 8.3 • 10^^ erg s~^ 

^ in the case of ionized hydrogen, deflection time for e~ — e~ collisions 
electron velocities 



and the energy was injected into two spherical regions of 
radius 13 kpc that lie at a distance of 30 kpc on either side 
of the centre of the cluster. This was done to simulate a 
bipolar outflow that is typically seen in clusters that show 
signs of AGN activity. The energy is injected by increasing 
the internal energy of the gas inside the injection region. 
We do not alter the velocity of the gas inside the injection 
region. The injected energy is uniformly distributed inside 
each injection region. The results are insensitive to details 
of the spatial profile of energy injection inside the injection 
regions. 

3.3. Dissipation of the injected energy 

The main transport processes responsible for the dis- 
sipation of the injected energy that we consider in this 
paper are viscosity and conductivity. Below we discuss 
their relative importance in the case of an unmagnetized 
plasma. 

Viscosity of the gas is due to off-diagonal terms in 
the stress tensor in the Navier-Stokes equations and, as 

such, is related to the momentum transport. /From or- 
der of magnitude estimates we get that viscous forces 
^ fiVgas/l'^; where Vgas is the typical gas velocity, I is the 
lengthscale over which the gas properties change and /x 
is the coefficient of viscosity. This coefficient is approx- 
imately given by /i ^ ViXiUnin, where n is the number 
density of gas particles, Vi is a typical velocity of species 
i (either electrons e or protons p) and \i and mi are their 
deflection length and mass, respectively. The deflection 
length is independent of mi^ (e.g., Frank, King and Raine 

1995). Since Vi am^ , this implies that the coefficient 

of viscosity is /i oc ^ rfhj'^ and that the transport of 
momentum and viscosity are dominated by protons. 

Both electrons and protons carry the same amount 
of energy. As electrons move faster than protons by 

{m.p/m,:!)^/'^ , the transfer of heat is mostly due to elec- 
trons. Thus, conductivity is dominated by electrons and 

the coefficient of conductivity is k oc rrie^^^. 

Let us now compare the relative contribution to dis- 
sipation from the above two processes. ^From order of 
magnitude estimates (c.f. equations 1 and 2) we have that 
Qcond ~ (K/fcB)eint/^^ and Qvisc ~ (M/mp)ekin/^^, where 
eint and ekin are the internal and kinetic energies of the 
gas and fc^ is the Boltzmann constant. Equipartition en- 
sures that the kinetic and internal energies are comparable. 
Thus, q = Qvisc/Qcond ~ {kB/mp){iJ,/K) oc {me/mpY/^. 
In an idealized case of plane linear waves in a uniform 
backround and without gravity the ratio of viscous to con- 
ductive dissipation turns out to be 



q = 



3 K 



(6) 



where Cy and Cp are the specific heat at constant volume 
and pressure, respectively (Landau & Lifshitz, 1997). Sub- 
stituting values appropriate for a fully ionized hydrogen 
plasma q = 4.2(me/mp)^/^ 0.1 (Braginskii 1958, Spitzer 

is shorter than for p — p collisions but this is compensated by higher 
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1962). Thus, in the average sense, the effect of conductiv- 
ity should dominate over viscosity in terms of the energy 
dissipation. In other words, this says that the ratio of the 
Peclet number to the Reynolds number is small. In the 
presence of magnetic field, both, viscosity and conductiv- 
ity are thought to be suppressed. However, the Larmor 
radius of ions is (mj/me)^/^ times larger that the Larmor 
radius of electrons of the same temperature. This can in- 
crease the Peclet number relative to the Reynolds number 
but their ratio may to exceed unity. 

We point out that, in this idealized case, the dissipa- 
tion length due to conductivity is shorter than that due 
to viscosity by a factor of q~^. However, wc also note 
that viscosity still plays a role in "diffusing" the gas mo- 
mentum by exerting viscous stress forces and affects the 
overall dynamics of the gas. Also, as demonstrated be- 
low, the regions where most of the conductive and viscous 
dissipation take place do not have to be spatially overlap- 
ping, i.e., some regions can be dominated by viscous or 
conductive dissipation. 

The dissipation of mechanical energy injected by the 
AGN can be estimated as follows. Assuming that the tem- 
perature varies smoothly with position in the fluid, the 
rate of change of mechanical energy can be written as 

+ I KVT-dS+ [ ^{VT)^dV (7) 
Js Jv T 

For random temperature fluctuations and |VT| ^ AT/L, 
where AT is the mean temperature change over some typ- 
ical length I/, the second term on the right hand side is 
small compared to the last one if we take the spatial aver- 
age of equation 7. In particular, for temperature fluctua- 
tions, such as waves, that occur in a constant background 
temperature, this term will vanish when averaged over one 
wavelength. In other words, if the temperature gradient 
is dominated by temperature changes on small scales, wc 
may locally approximate the dissipation of mechanical en- 
ergy by the last term in equation 7. This is the approach 
that we adopted when comparing the dissipation due to 
viscosity and conductivity with the radiative cooling. 

The simulations of the effects of viscosity and conduc- 
tivity in three dimensions require high spatial resolution. 
This can be achieved at the expense of a relatively short 
simulation time (a few full AGN activity cycles). We 
point out that the timestep imposed by the transport co- 
efficients on the simulation scales as (A.t)^/(k, /i), where 
Ax is the simulation resolution. As this timestep scales 
more strongly with Ax than the standard hydrodynam- 
ical Courant condition, and because the transport coef- 
flcients depend strongly on temperature, the constraints 
on the timestep are more stringent that the ones obtained 
from the standard condition. This is why we chose to re- 
sort to comparing instantaneous rates of heating instead 
of evolving the system for a longer time at lower spatial 
resolution as the latter could prevent us from capturing 
essential physics. However, we note that the actual evo- 
lution of the gas in the simulations, including the energy 
transfer and dissipation, does not use the above approxi- 



mation. We also point out that the actual heating of the 
cluster core is greater than that estimated from this pre- 
scription. This is because the spatially and temporally 
averaged temperature gradient in a cooling flow cluster is 
positive within the cooling radius and the heat transfer 
from the hot outer layers will take place. Thus, the es- 
timated conductive dissipation can be considered to be a 
lower limit on the actual dissipation due to conductivity. 

4. SYNTHETIC CHANDRA OBSERVATIONS 

We have performed synthetic X-ray observations of the 
grid at two stages of bubble evolution. These observa- 
tions simulate a 200ks ACIS-S3 exposure of the (800A;pc)^ 
region centered on the bubble origin. The images are gen- 
erated by first calculating a MEKAL (Mewe et al. 1985) 
cmissivity at each grid point in three energy bands, then 
integrating along each linc-of-sight through the simulation 
box in the optically thin limit using the SYNTH code (e.g., 
Tregillis et al. 2004). The resulting X-ray surface bright- 
ness image in each band is modified by galactic absorption, 
adjusted for the assumed redshift of the simulated cluster 
(z=0.0183, Perseus) then processed through the Chandra 
ACIS-S3 response function to generate instrument count 
rates at each pixel of the image. The effect of vignetting is 
included and corrected for in the final images. In this case 
the image plane is larger than the dimensions of ACIS-S3, 
so we tile the chip accordingly to cover the image plane. 
Since the angle subtended by the smallest voxel in the sim- 
ulation when placed at Perseus' redshift is approximately 
7.5", no modeling of the PSF is necessary. 

To simulate the effect of the instrumental and X-ray 
background, we have added the ACIS-S D period sky back- 
ground to the resulting images. The background event 
file is sorted to generate images in each of the three X- 
ray bands, then binned to match image resolution, and 
matched to the exposure of the synthetic observations. To 
simulate the effect of background subtraction, we have 
simply mirrored each of the three background images 
about the vertical chip axis and subtracted the mirrored 
image in each case. The net counts added and subtracted 
to the image then are the same, but the counts in each 
pixel arc not identical, resulting in the appearance of a 
background-subtracted image. The final images have the 
visual appearance of a Chandra X-ray image. More im- 
portantly, in this way we can quantify the expected error 
level associated with the features in the images. 

The final images are equivalent to exposure corrected, 
background subtracted ACIS-S3 images of the simulated 
volume. In addition, for the images in each of the three 
energy bands we have performed an unsharp masking to 
generate an additional image. This image results from 
subtracting an image smoothed with a 4 pixel gaussian 
kernel from the original image. This procedure enhances 
the appearance of fluctuations about the mean. 

5. RESULTS 

Figure 3 shows three slices through the cluster center 
displaying the gas temperature at the time of 140 Myrs 
(file 0142) after the start of AGN activity. The panels 
have the size of 2.8 Mpc, 800 kpc and 300 kpc on a side 
from left to right, respectively. At this stage in the evolu- 
tion one can clearly identify these bubbles that have not 
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Fig. 3. — Temperature distribution in a slice through the cluster center. The box sizes in the panels are 2.8 Mpc, 800 kpc, and 300 kpc 
from left to right, respectively. All panels correspond to the time of 140 Myr after the initial AGN outburst. Both young bubbles still in the 
inflation phase as well as older buoyantly rising bubbles can be seen. Note nearly spherical weak shocks surrounding the bubbles. All plots 
are logarithmic and the color bars show the logarithm of temperature in K. 



yet reached pressure equilibrium with their surroundings 
and that are still expanding nearly spherically into the 
ambient medium (bubbles close to the injection region) as 
well as older bubbles that have evolved into mushroom- 
type clouds. As a result of the rapid inflation, the younger 
bubbles produce a weak shock wave that can be seen to 
travel outward. Figure 1 also shows that the cluster is 
quite dynamic, as it shows significant substructure such 
as clumps and shock fronts. 

Figure 4 shows the density distribution in the cluster. 
This figure corresponds to the same time and box sizes 
as Figure 1. As in Figure 1 one can clearly identify the 
young bubbles that are still in the expansion phase and 
the older, buoyant bubbles. As expected in the case of 
weak shocks, the waves are weaker in the density maps 
than in temperature maps but are still visible. 

Although the detached bubbles shown in Figures 3 and 
4 do not seem to be disrupted (their detailed appear- 
ance depends on the exact position of the cross-section 
plane), they eventually do tend to mix with the intraclus- 
ter medium late in the bubble evolution. Observations 
indicate that the bubbles maintain their identity as they 
move away from cluster centers. Most observations detect 
bubbles in the early stages of their evolution where vari- 
ous instabilities may not have had enough time to disrupt 
them. However, if the real bubbles indeed preserve their 
identity at larger distances then this effect is not explained 
by our simulations. This discrepancy could then be ex- 
plained if either the bubbles themselves were viscous (as 
discussed by Reynolds et al. 2004) or if they were magne- 
tized in which case Rayleigh- Taylor and Kelvin-Helmholtz 
instabilities on the bubble surfaces would be suppressed. 
We note that even initially weak magnetic fields (i.e., dy- 
namically unimportant ones) could be amplified around 
the bubbles and effectively prevent instabilities from hap- 
pening (Jones & DeYoung 2005). 

Figure 5 presents the dissipation rates corresponding to 
the time of 80 Myr after the initial AGN outburst. The 
left panel shows the ratio of the viscous dissipation to the 
radiative cooling. As can clearly be seen, the weak shocks 



present in the density and temperature maps are the sites 
of enhanced dissipation. Some of the strongest waves have 
dissipation rates comparable to the radiative cooling rates. 
The middle panel shows the ratio of the conductive dissi- 
pation and the cooling rate. A comparison of this and 
the left panel shows that the conductive dissipation in the 
shocks on average appears to be higher. This is also visible 
in the right panel that shows the ratio of the conductive 
and viscous dissipation rates and especially in the profiles 
presented below. 

Figure 6 is analogous to Figure 5 but corresponds to a 
later time (140 Myr) after the onset of AGN activity. The 
dissipation patterns presented in this figure correspond to 
the same time as the temperature and density maps in 
Figure 1 and 4, respectively. Comparison of Figures 5 and 
6 shows that the dissipation patterns moved away from 
the center. The typical speed of these patterns is of or- 
der of the sound speed in this cluster. This reinforces the 
interpretation of these features as sound waves or weak 
shocks. 

Figures 7 and 8 present the profiles of viscous and con- 
ductive dissipation for the time of 80 Myr and 140 Myr 
after the first outburst, respectively. The profiles were ex- 
tracted along the horizontal lines intersecting the centers 
of the dissipation maps shown in Figures 5 and 6. Left 
panels present the profiles of the ratios of viscous dissipa- 
tion and cooling rate. The middle ones are for conductive 
dissipation and the right panels show the profiles of the 
ratio of conductive and viscous dissipation. It is interest- 
ing to note that the conductive-to-viscous ratios seem to 
be 10 which is consistent with simple analytic argument 
presented in Section 3.3. 

Figures 7 and 8 indicate that there is a large degree of 
scatter in the ratio of the heating rates. The deviations 
of the simulated ratios from the estimated constant ratio 
given by Equation 6 are mostly due to the fact that the 
density and temperature fluctuations are nonlinear (the 
gas may be effectively "swept" as a result of energy injec- 
tion) and due to the the fact that bulk gas motions are 
present in the cluster. We note that the offsetting of the 
effects of cooling is helped by the fact that bubbles tend 
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Fig. 4. — Same as Fig. 3 but for physical gas density. 




Fig. 5. — Dissipation patterns. All panels are 300 kpc a side and correspond to 80 Myr after the initial outburst of the AGN. Left panel 
shows the logarithmic ratio of the viscous dissipation and the radiative cooling rate. Middle panel presents the logarithm of the ratio of 
the conductive dissipation to the radiative cooling rate and the right one presents the logarithm of the ratio of the conductive and viscous 
dissipation rates. 




Fig. 6. — Same as Figure 5 but for the simulation time of 140 Myr. 
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Fig. 7. — Profiles along the horizontal lines through the center of the dissipation patterns in Figure 5. 




Fig. 8. — Profiles along the horizontal lines through the center of the dissipation patterns in Figure 6. 



to transfer momentum to the ambient gas and disperse it. 
The gas that gets elevated and disperssed can be heated 
more efHcicntly. 

Figures 9 and 10 show X-ray images of the chistcr heatcnl 
by AGN at the time of 80 Myr and 140 Myr, respectively. 
The first and the second row show the maps of X-ray sur- 
face brighncss and unsharpcd X-ray images, respectively. 
The columns are for [3.5-7.0] keV [1.5-3.5] keV and [0.3-1.5] 
keV from left to right. These figures show clear evidence 
for X-ray featureas at the spatial locations of the waves 
resulting from expansion of the bubbles. The bubbles and 
waves can be clearly seen in the earlier stages in the AGN 
evolution. All features are more easily detectable in the 
softer X-ray bands. 

Figure 11 shows a closeup of the unsharp masked image 
in the soft band corresponding to the time of 80 Myr (bot- 
tom right panel in Figure 9). Apart from the waves one 
can easily identify a bullet-like feature wrapping around 
one side of the maximum in density distribution suggest- 
ing a motion to the right. 

6. SUMMARY 

We have simulated the effects of heating by active galac- 
tic nuclei on the intracluster medium including, both, the 
effect of viscosity and conductivity. This has been done by 
a high resolution 3D AMR simulation of a galaxy cluster 
that has been extracted from full a cosmological SPH sim- 
ulation. Our simulations demonstrate that conductivity is 
likely to play an important role in dissipating mechanical 
energy injected by the AGN. We have shown that if the 
suppression of conductivity and viscosity is of comparable 
magnitude, then the effect of conductivity on dissipation is 
likely to exceed that due to viscosity. Nevertheless, viscos- 
ity may play a significant role in affecting the dynamical 
evolution of the buoyant bubbles especially if suppression 



of conductivity is higher than that of viscosity. The precise 
values of suppression factors are unknown as yet but we 
hope that future observations of the ICM would be able 
to put c:onstraints on their values. Our simulations were 
performed for a moderate mass cluster. As higher mass 
clusters have higher temperatures and densities and the 
opposite is true for lower mass ones, we expect that quali- 
tatively the trends found in our simulations could hold for 
clusters different than the one considered here. However, 
we note that the details of feedback are sensitive to the 
amount of the injected energy that may vary depending 
on the cluster mass (Roychowdhury ct al. 2004). They 
also depend sensitively on the exact values of density and 
temperature as they enter in high powers in the expres- 
sions for cooling and dissipation rates. We also discuss 
the detectability of the X-ray features generated by AGN 
outbursts. Our synthetic data takes into account all in- 
strumental features of the Chandra Observatory and all 
sources of noise such as X-ray background. We show that 
bubbles inflated by AGN and the waves generated by the 
expansion of the bubbles can be seen in these synthetic ob- 
servations. The detactability of tliese features depends on 
the stage in the AGN evolution (with bubbles more easily 
detectable in the earlier stages) and the Chandra energy 
band observed (with stronger features in softer bands). 
Other features such as bullets due to substructure mo- 
tions are also detectable in the simulated X-ray emissivity 
maps. These are similar to that observed in IE 0657-56 
(Markevitch et al. 2002) 
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Fig. 10. — Same as Figure 9 but for the time of 140 Myr. 




Fig. 11. — Closeup of the bottom right panel in Fig. 9. Both the waves and the "bullet-like" feature are visible. 
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